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Abstract 

We examine whether the cosmologies with varying speed of light (VSL) are com- 
patible with the second law of thermodynamics. We find that the VSL cosmology with 
varying fundamental constant is severely constrained by the second law of thermody- 
namics, whereas the bimetric cosmological models are less constrained. 
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1 Introduction 



Variable-Speed-of-Light (VSL) cosmological models were proposed |], @] as an alter- 
native to inflation || |5| for solving the initial value problems in the standard Big 
Bang model. It is assumed in the VSL models that the speed of light initially took 
a larger value and then decreased to the present day value during an early period of 
cosmic evolution. VSL models have attracted some attention, because not only various 
cosmological problems that are solved by the inflationary models but also the cosmo- 
logical constant problem can be solved Q, |, g % |, |, |, |, § § by VSL 
models. Furthermore, the recent study of quasar absorption line spectra in comparison 
with laboratory spectra shows that the fine structure constant a = e 2 / (AnTic) varies 



over cosmological time scales |Tq, 17, 18 1, indicating that the speed of light may indeed 



vary with time. Also, it has been shown |19], |20, EH], E2l [23], E3, EH |26| that brane world 



models, which have been in vogue recently, manifest the Lorentz violation, which is 
a necessary requirement for the VSL models. (Cf. The recent work |27J studies the 



experimental limits which are permitted for the graviton's speed in the brane world 
scenarios.) 

It is the purpose of this paper to examine the compatibility of the VSL models with 
the second law of thermodynamics. (Cf. The previous related work can be found in Ref. 
||28|| .) Recently, there has been active interest in holographic principle in cosmology, 
after the initial work by Fischler and Susskind (FS) f29|. The cosmological holographic 
bound originally formulated by FS had a problem of being violated by the closed 
Friedmann-Robertson- Walker (FRW) universe. Later works attempted to circumvent 
such a problem through various modifications. In particular, it was proposed in Refs. 
fl30| , PH] that the FS holographic bound has to be replaced by the generalized second 
law of thermodynamics. The generalized second law states that the total entropy S of 
the universe should not decrease with time during the cosmological evolution: dS > 0. 
In order to be compatible with the holographic principle, the VSL cosmological models 
therefore have to obey the generalized second law of thermodynamics. In section 2, 
we consider the original VSL model, where the speed of light c in the action is just 
assumed to vary with time. In section 3, we consider the bimetric cosmology of Clayton 
and Moffat. 



2 VSL Cosmology with Varying Fundamental Con- 
stant 

First, we consider the original VSL cosmology [0, 0, in which a fundamental constant 
c of the nature is just assumed to vary with time during the early period of the cosmic 
evolution and thereby the Lorentz symmetry is explicitly broken. In such VSL theories, 
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it is postulated that there exists a preferred Lorentz frame in which laws of physics 
simplify with the action taking a standard form with a constant c replaced by a field 
c(x /1 ), the principle of minimal coupling. Namely, the action in the preferred frame 
takes the form 



S = aTx 



c A 



(1) 



where C c controls the dynamics of c and C is the action for the fields in the universe. 
It is required that C c should be explicitly independent of the other fields, including 
metric, so that the principle of minimal coupling continues to hold for the equations 
of motion. 

The general metric ansatz for a 4-dimensional homogeneous and isotropic universe 
is given by the following Robert son- Walker metric: 

g^dx^dx 11 = —c 2 dt 2 + a 2 ^ijdx % dx\ (2) 

with the time- varying c. Here, a(t) is the cosmic scale factor and r jij(x k ) is given by 

/ \ —2 dv 2 

li3 dx l dx J = (l + \ 6 mn x m x n ) 5 ij dx t dx ] = 1 _ kr2 + r 2 dfl 2 2 , (3) 

where k = —1, 0, 1 respectively for the open, flat and closed universes. 

With an assumption of the principle of minimal coupling, the Einstein equations 
with the metric ansatz (Q) lead to the Friedmann equations: 

d\ 2 8nG c 2 . kc 2 
ay 6 6 cr 

where the overdot denotes derivative w.r.t. t. From the Friedmann equations, we 
obtain the following generalized conservation equation: 

/ p \ a 3kcc cc , 
p + 3[p+ ^) - = —— - -—A. 6 
V c 2 J a AirGa 2 AnG 

We now study thermodynamics of the VSL cosmology. We assume that the uni- 
verse satisfies the first law of thermodynamics. When applied to the comoving volume 
element of unit coordinate volume and physical volume v = a 3 , the first law of ther- 
modynamics takes the form: 

Tds = d(pc 2 v) + pdv, (7) 
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where s = s(v,T) is the entropy density of the universe at temperature T within the 
volume v, and p = p{T) and p = p(T) are the mass density and the pressure of matter 
in the universe. In this paper, we assume that c is a function of a, just as in Ref. ||. 
Since v = a 3 , we have c'(a) = 3a 2 dc/dv = 3v 2 ^ 3 dc/dv, where prime denotes derivative 
w.r.t. a. Then, Eq. (|7|) can be rewritten as 

Tds = c 2 vdp + ^pc 2 + p + 2pcv — ^ dv 

= c 2 vdp + (pc 2 + p + \pcdv 1/3 ) dv. (8) 

So, partial derivatives of s(v, T) are given by 

9S{V ' T) =^(pc 2 +p+lpcc'v^ 



dv 



ds(v, T) c 2 v dp 
dT = ~YdT' 

From the integrability condition d 2 s/(dvdT) = d 2 s/(dTdv), we obtain 

^ (pc 2 + p+ Ipcc'v 1 /^ = ^— [(pc 2 + p)v . 
Making use of this equation, we can put Eq. ([/J) into the following form: 

ds = d 



dp 
dT 



(9) 
(10) 

(11) 



2v 2 dc lrT1 



(12) 



from which we see that the usual Euler's relation s = ^(pc 2 + p) does not hold for the 
VSL theories with varying fundamental constant. 

To obtain the time derivative of s, we express the conservation equation (|(J) into 
the following form, making use of Eq. (|TT|) : 



d 
dl 



T 



(pc 2 + p) 



c cc 3kcc 



c 2 v 



dc o T 



+ 2pc—v 



dv T 2 ' 



(13) 



Eq. (|T2D along with Eq. (|T3|) yields Q 



2p 



-A 



3kc 2 



cca 



AnG AirGa 2 



(14) 



This equation would also have been obtained directly from Eqs. 



and (0). 



2 The previous related work |28| does not take into account the effects on s of the modification of 
the usual Euler's relation due to the time varying c and the nonzero cosmological constant A. 
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Since c < for the VSL cosmology with varying fundamental constant, the terms 
in the parenthesis of Eq. (|14]) have to be non-positive in order to be compatible with 
the second law of thermodynamics dS > 0: 

2p-— A + — — <0. 15 
AttG AnGa 2 

When the cosmological constant is non-positive (A < 0), this condition can never be 
satisfied by the flat (k = 0) and the closed (k = 1) universes. Although it can be 
satisfied by the open universe (k = —1), the condition is very restrictive about the 
possible type of matter in the universe and the time variation of c. On the other 
hand, when the cosmological constant is positive (A > 0) and is large enough, the 
restriction becomes less severe. We now examine the explicit condition under which 
the constraint ([15]) can be satisfied. Assuming that c varies gradually like c(t) = CqoJ 1 
and the equation of state of the form p = (7 — l)pc 2 , by solving Eq. ([]) one obtains 

B 3kc 2 na 2( - n -^ Anc 2 a 2n 

P ~~^ + 4fcG(3 7 + 2n - 2) ~ 4ttG(37 + 2n) ' ^ 

where B is a positive constant. So, the constraint fljjD reduces to 

2B , 3(3 1 + An -2)kc 2 < (3 1 + An)c 2 A 



a 3 7 +2n 4 vr c(3 7 4. 2n - 2)a 2 ~ 4vrG(37 + 2n) 

First of all, we see from this that for the flat universe (k = 0) the constraint is always 
violated for some values of a. [When 37 + 2n = 0, the last term on the RHS of Eq. 
(16f) behaves with a as ~ a -37 In a and therefore the RHS of Eq. ( |17D behaves with a 
as ~ In a.] In order for Eq. QT7D to be satisfied for any values of a, the LHS has to 
have the maximum value and its value has to be not greater than the RHS. From such 
condition, we obtain the following constraint on the constants n and 7: 

A;(3 7 + An - 2) > 0, (3 7 + 2n)(3 7 + 2n - 2) < 0, (18) 

along with the following lower limit on the cosmological constant: 

3(3-i + An - 2)kc 2 Q 



A> 3(3 7 + 4n-2)fe 
3 7 + An 



4tt(3 7 + 2n)(3 7 + 2n - 2)GB 



3t+2u-2 



(19) 



One can always choose the integration constant B in such a way that the bound ( JL9D 
is compatible with the observed value of A. However, the constraint fll8|) severely 
restricts the allowed value of 7. For the closed universe (k = 1), Eq. ( |i~8| ) implies 
7 < 2/3 and therefore the radiation-dominated universe (7 = 4/3) and the matter- 
dominated universe (7 = 1) are not allowed. As for the open universe (k = —1), Eq. 
([[ID leads to the less severe constraint 7 > —2/3, which allows the 7 = 4/3, 1 cases. 
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3 Scalar- Tensor Bimetric Cosmology 



In this section, we consider the bimetric VSL model, proposed by Clayton and Moffat 
TIJ. The bimetric models achieve time- variable speed of light in a diffeomorphism 



invariant manner and without explicitly breaking the Lorentz symmetry by introducing 
bimetric structure into spacetime. (Cf. It is recently found out in Ref. [02] that the 



fine-structure constant a = e 2 / (4irhc) in the bimetric models is constant in spacetime 
although the speed of light varies with time, due to the compensating time- variation 
of the electric charge.) It is usually assumed in the bimetric models that graviton and 
the biscalar (or the bivector) propagate on the geometry described by the "gravity 
metric", whereas all the matter fields (including photons) propagate on the geometry 
described by the "matter metric" . In the case of the scalar-tensor bimetric model, the 
gravity metric and the matter metric g^ are related by the biscalar field $ as 

g^u = g^u - Bdfidv®, (20) 

where a dimensionless constant B is assumed to be positive. Since these two metrics 
are nonconformally related, a photon and a graviton propagate at different speeds. The 
action for the scalar-tensor bimetric model has the form 

S = J d 4 xV=g~ (K -2A) + U + J d 4 x^g~£ mat , (21) 

where £ ma t is the Lagrangian density for matter fields and the Lagrangian density 
for the biscalar is given by 



C* = ~g' w d li $d v Q-V{$). (22) 



The gravity metric for the universe has the form 

g tw dx p 'dx v = —c 2 dt 2 + a 2n fijdx % dx^ , (23) 

with constant speed of graviton c grav = c and 7^ given by Eq. ([3]). Due to the require- 
ment of homogeneity and isotropy of the universe, the biscalar field $ is independent 
of the spatial coordinates x l . According to Eq. (|20|), the matter metric is therefore 
given by 

g^dx^dx" = -(c 2 + B&)dt 2 + a 2 'y ij dx i dx i , (24) 
where the overdot stands for derivative w.r.t. t. So, the speed of photon c p h = 



cy 1 + B& 2 /c 2 = c\fl varies with t, taking larger value than c grav = c while $ 7^ 0. 

In obtaining the energy-momentum tensor for the purpose of deriving the Einstein's 
equations, one has to keep in mind that matter fields and biscalar are coupled to 
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different metrics. Since g^ v is the physical metric for the matter fields, the energy- 
momentum tensor for the matter fields are defined in terms of g^: 

fiXV ^ 2 5(V^gC mat ) = ^ + ^ ^ + ^ (25) 



89, 

where p and p are the mass density and the pressure of the matter fields and is the 
four- velocity of matter perfect fluid normalized as g flu U fl U L ' — — 1. Since the nonzero 
component of the four- velocity vector in the comoving coordinates is U* = l/c p h, the 
nonzero components of the energy-momentum tensor for the matter fields are 

ftt = p fij = V / 2 gx 

er 

On the other hand, since the biscalar field is coupled to the gravity metric g^, its 
energy- momentum tensor is defined in terms of g^ v : 



7T = 2_ S (V~9£*) = gtWgvPd^dp® - \g^d a §d a <$> - V{<S>)g^ 
V 9 v9nv 

= (p$c 2 +p„)u»U»+p^ u , (27) 



where the four- velocity for the biscalar is normalized as g llv U^U l/ = 1, so its nonzero 
component is U* = 1/c. The mass density and the pressure of the biscalar field are 
therefore 

1$2 \ 1 i$2 



p*=[«^r + v)-z, p* = ~-j-v. (28) 



Taking the variation of the action S w.r.t. the metric, we obtain the Einstein's 
equations 

Q»v + A g ^ = (29) 

c 4 

where Q^ v is the Einstein tensor for g^ and the energy-momentum tensor has the 
form 

= rr + T^^2. (30) 
The nonzero components of are 

r« = P$ + p^, = f P * + p^) V- (si) 

The Einstein's equations lead to the Friedmann equations 

a\ 2 8ttG c 2 . kc 2 , nn . 
- =—p + -A~— , (32) 
aj 6 6 a z 
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a = -^l P + 3 ?J + 3 A ' (33) 

where 

p = p$+p— , p = p $ +p^. (34) 

c c 

From these Friedmann equations, we obtain the following conservation equation 

p + 3fp+4)- = 0. (35) 
V c z J a 

The factors of c p h/c grav = c p h/c in Eq. ( P^ ) can be understood from the fact that 
definition of the energy-momentum tensor depends on the choice of metric. Generally, 
the following two energy-momentum tensors, associated with the same Lagrangian 
density L but defined w.r.t. the two different metrics g^ and g^, 

T» v = 2 6L = 2 6L (*R\ 



are related to each other as 



i grav _ / — -"'ph _ J ph- l^'J 



"9 c g 

Here, the subscripts 'grav' and 'ph' signify that the quantity under consideration is 
defined w.r.t. the gravity and the matter metrics, respectively. The mass densities and 
the pressures in the two different definitions are defined by 

-^grav = (PgravCg rav + Pgrav) U^U + Pgravfi 1 ^ > 

3ff = (pphc; h +p ph )f>^+p ph ^, (38) 

where the four-velocities are normalized as g^uU^U" = —1 and g tlv U ll tj v = —1. From 
Eqs. (|3~7p and (|38|), we see that the mass densities and the pressures in the two different 
definitions are related as 

Pgrav Pph) Pgrav Pph- (39) 

Cgrav Cgrav 

The factor of c grav /c p h, multiplying the matter fields mass density and pressure in Eq. 
(|54|), arose due to the fact that the matter fields mass density and pressure, which 
are defined w.r.t. the matter metric, has to be transformed to the 'gravity metric' 
quantities, since the Friedmann equations are defined w.r.t. the gravity metric. (Cf. 
In the Einstein's equations (|29|), the total energy- momentum tensor T M1/ is defined 
w.r.t. the gravity metric, i.e., T^ v = where Lagrangian density L is the 

sum of the matter fields Lagrangian density L mat = \/—gC mat and the biscalar field 
Lagrangian density L$ = ^—gC§.) 
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We now study thermodynamics of the bimetric VSL cosmology. Since we are con- 
sidering the frame associated with the gravity metric of the form (|23|) , the physical 
quantities in the thermodynamic laws should be 'gravity metric' quantities. The first 
law of thermodynamics, applied to the comoving volume element of unit coordinate 
volume and physical volume v = a 3 , therefore takes the form 

Tds = d{pc 2 v) + pdv, (40) 



where p and p are given by Eq. (|34j) . Note, the square of the speed of graviton c grav = c 
multiplies p to yield the energy density, because we are considering the frame associated 
with the gravity metric. From Eqs. (|35|) and ([40"D , we obtain ds/dt = 0, namely the 
total entropy S = s J dx 3 ^/—^ of the universe remains constant during the cosmic 
evolution. Therefore, unlike the VSL cosmology with varying fundamental constant, 
considered in the previous section, the second law of thermodynamics dS > is always 
obeyed regardless of values of k. 

We comment on thermodynamics associated with the matter fields, only. Since it is 
assumed that the matter field action is constructed out of g^, the equations of motion 
of the matter fields imply the conservation law for the matter fields energy-momentum 



tensor 12 : 



V M T^ = 0, (41) 

where V M is the covariant derivative defined w.r.t. g^ u . The t-component V M 2j = 
of the conservation equation, where = T w g pv , takes the following form: 

* + 3 (* + *-)° = -2^. (42) 

From this equation we see that matter fields are created while the speed of photon 
decreases with time to the present day value c. In the frame associated with the 
matter metric, the first law of thermodynamics takes the form: 

Tds mat = d(pc 2 vh v) + pdv, (43) 

where s ma t is the entropy density associated with the matter fields. From Eqs. (|42"1) 
and (^), we see that s mat = 0, namely that the entropy for the matter fields remain 
constant despite that the matter fields are created while c p h < 0. This apparent 
paradox can be understood from the fact that in the matter metric frame the past 
light cone contracts and thereby less information is collected by the observer. Next, we 
consider the frame associated with the gravity metric. In this frame, the mass density 
and the pressure of the matter fields are transformed to pc^/c and pc p h/c, and the 
square of the speed of graviton c should be multiplied to the mass density for obtaining 
the energy density. So, the first law of thermodynamics takes the form: 

Tds m3ut = d(pcc ph v) + ^-pdv. (44) 
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From Eqs. and (44), we have 

_ ( (? ph -<? a _ c ph \ cc ph a 3 

So, in the gravity frame, the entropy for the matter fields varies with time while 
Cmat 7^ 0. This time variation of the entropy of the matter fields is due to the exchange 
of entropy with the biscalar sector, since we have seen that the total entropy density s 
remains constant in the gravity metric frame. 

In the above, we considered the equations in the comoving frame for the gravity 
metric g^ v . Since all the matter fields in the universe are coupled to the matter metric 
g^ u , it would be more natural to consider the comoving frame for the matter metric in 
studying the expansion of the universe. By defining the cosmic time r for the matter 
metric in the following way 

dr 2 = (l + B<b 2 /c 2 )dt 2 } (46) 

we can bring the matter metric into the following standard comoving frame form for 
the Robert son- Walker metric: 

g^dx^dx" = -c 2 dr 2 + a 2 (T)-f ij dx i dx j . (47) 



In this new frame, the gravity metric (|23|) takes the form 



g^dx»dx u = -(c 2 - B<$> 2 )dr 2 + ^(r^dx'dx 1 , (48) 

where the overdot from now on stands for derivative w.r.t. r. So, in this new frame, 
a photon travels with a constant speed c p h = c and a graviton travels with a time- 



variable speed c grav = y c 2 — B<& 2 = c/\/l, taking smaller value than c while $ ^ 0. 
Note, / = 1/(1 — B$ 2 /c 2 ) when the overdot stands for derivative w.r.t. r. 

In obtaining the Friedmann equations in the new frame, we do not just apply 
the change of time coordinate JIB] ) in the Friedmann equations ( ft^ ) and ([£J) in the 
old frame, unlike the previous works on bimetric VSL cosmology. The reason is that 
the definitions for the mass density and the pressure depend on the choice of time 
coordinate. We consider the Einstein's equations ( |29"D with Q^ u now being the Einstein 
tensor for the gravity metric given by Eq. ©. The energy-momentum tensor is 
still given by Eq. ([30]) but now with 

f" v = (pc 2 + p) U»U V + pg^, (49) 

n u = ( P$ 4 av + p*) uw + p <r , (so) 
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where the nonzero components of the four- velocities are U* = 1/c and U* = 1/cgrav 
Note, the mass densities and the pressures in Eqs. (|49| ) and fl50"l) are different from 
those in Eqs. fl25|) and (0), although they are denoted with the same notations. In 
the comoving frame for the matter metric, the Friedmann equations therefore take the 
forms 



a) 3 3 a 2 



a Cgrav CL AnG ( , P \ , C f 

O Cgrav Q ^ V ^grav 



P + 3^- + -^A, (52) 



where 

c c 
p = p$ + p , p = p^+p . (53) 

Cgrav Cgrav 

Note, although we are now considering the comoving frame for the matter metric, the 
total energy momentum tensor is still defined w.r.t. the gravity metric. It is just 
that the time coordinate r in the gravity metric is the comoving frame time coordinate 
for the matter metric. So, mass density and pressure of matter fields still have the 
factor of c/c g rav in Eq. (fxf). From these Friedmann equations, we obtain the following 
conservation equation 

. . / P \ & 3 Cgrav f Q'\ CgravCgrav . 3 &Cg r avCgrav fr a\ 

p + 3[p+— — I - = -r~^ — \ % A + ~a — — 9"^ — ■ ( 54 ) 



^grav / 



We now study the compatibility of the bimetric VSL cosmology with the second 
law of thermodynamics. Although we are now considering the comoving frame for the 
matter metric, we first study thermodynamics in the frame of the gravity metric, since 
the conservation equation (54) is expressed in the gravity metric frame. The first law 



of thermodynamics, applied to the comoving volume element of unit coordinate volume 
and physical volume v = a 3 , takes the form 

Tds = d(pCg rscv v) + pdv. (55) 

From Eqs. ([>|) and (|55D, we obtain 

TS = 4pa 3 CgravC g rav, (56) 



where we made use of Eq. (plf ) to simplify the RHS. So, the total entropy S = 
s J dx 3 ^—"f increases with time while c gra v — c/ \/T increases to the present day value. 
The second law of thermodynamics is therefore always satisfied for any values of k. 
Next, we consider the second law of thermodynamics in the frame of the matter metric. 
The mass density and the pressure in the matter metric frame are given by pc grav /c 
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and pc grav /c, and the energy density is obtained by multiplying the mass density by c 2 . 
So, the first law of thermodynamics takes the form: 



Tds = d(pcc graiV v) + grav pdv. (57) 



This along with Eq. (|54|) leads to 

a = 3 ( °% av ~ ° 2 p- + p^-) CC ^ a \ (58) 

\ C C grav a Cgrav/ T 

So, the terms in the parenthesis of this equation has to be non-negative in order for 
the the condition for the second law of thermodynamics to be satisfied. Making use of 
the explicit expressions for c grav , p and p, one can put the condition in the form: 



1 I 



2 c 



2 



^ 2 -v + pVi - + ^4^ 2 + v + pc 2 Vi 



a 



1 I 



2 c 



2 



< 0. (59) 



Making use of the equation of motion for the biscalar field it can be shown that 
regardless of value of k this condition can be satisfied by the perfect fluid obeying the 
equation of state p = (7 — l)pc 2 with 7 < 2, provided the biscalar potential satisfies 
d T V(<$>) > 0. 

We discuss thermodynamics associated with matter fields, only. As before, the 
equations of motion for the matter fields imply the conservation law ( f4ip for the matter 
fields energy-momentum tensor. Since the speed of a photon takes a constant value 
c p h — c in the comoving frame for the matter metric, the t-component of the matter 
fields conservation equation takes the form: 

'P + 3 {P+^)1 = °- ( 6 °) 

Unlike the case of the comoving frame for the gravity metric, the matter fields are 
observed to be conserved. Eq. (|60| ) along with the first law of thermodynamics for the 
matter fields 

Tds mat = d(pc 2 v) + pdv, (61) 

yields s ma t = 0. So, time variation of entropy density s in the matter metric frame, 
as expressed in Eq. (|58|) , is all due to time variation of entropy density of the biscalar 
field, and there is no entropy exchange between matter and the biscalar field sectors. 

4 Conclusions 

We studied the compatibility of the VSL cosmological models with the second law of 
thermodynamics. We find that generally the VSL model with varying fundamental 
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constant is severely constrained by the second law of thermodynamics. For the model 
with the speed of light gradually varying as a power law of cosmic scale factor, only 
open universe with positive cosmological constant is allowed to contain cosmic perfect 
fluid satisfying reasonable equation of state such as dust and radiation. On the other 
hand, the bimetric cosmological models are less constrained by the second law of ther- 
modynamics. In the comoving frame associated with the gravity metric, there is no 
constraint on the model by the second law of thermodynamics. In the comoving frame 
associated with the matter metric, provided the biscalar potential satisfies a certain 
condition, regardless of the value of k, the second law of thermodynamics can be sat- 
isfied by perfect fluid satisfying the equation of state p = (7 — l)pc 2 with 7 < 2, which 
includes dust and radiation. 
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